TH associação dos estudantes do instituto superior técnico 
dezembro 1971 revista de engenharia 


VARÕES DE AÇO DE ALTA RESISTÊNCIA 


PARA BETÃO ARMADO TIPO A. 40 
(«Heliaço» e «SNT 40») 


TENSÕES DE SEGURANÇA 
Ra = 2400 kg/cm? 
e 
Ra = 4000 kg/cm? (SISMOS) 


Literatura e Tabelas de Cálculo 
fornecidas gratuitamente 


HELIAÇO, LDA. 


R. DE S. JULIÃO, 41, 4.º D. — LISBOA 


Tel.: 36 7817-3237 06-369477 


materiais 
novobra 


LISBOA - LEIRIA - LAGOA - GUARDA 


TUBOS DE BETÃO 


| A 
NÃ 
++ 


0 ARMADOS 

DE FABRICO INTEIRAMENTE AUTOMÁTICO 
STOCK PERMANENTE 

LISBOA -5 v 


LEIRIA Ponte « 
LAGOA Estrada do Carmo -Te 4 


GUARDA  Arritana - Tel 817 


MEXIA HEITOR E BRASÃO FARINHA 


aa 


a ASSOC dação dos estudantes do ins 


ISSO A 


LIMENTO Seb 


(22 


Fábrica: 
Outão — Setúbal 


Telef. 228824/5 


Escritório: 


KR, do Comércio, 56-3.º 


LISBOA 2 
Telef. 328201/2/3 


DEPARTAMENTO DE 


OT 
So 
a 


DO 


ENGENHARIA 
MECANICA 


— o — 


p= 


NUMERO 413 


DIRECTOR 
Guilherme Silva Arroz 


ADMINISTRADOR 
José Toscano 


CHEFE DE REDACÇÃO 
Francisco Pedroso 


SECRETÁRIO 
Josó Sá da Costa 


REDACTORES 


António Mouraz Miranda, António 
Vieira Pita, José de Oliveira Paulo, 
José Xambre Bento, Maria Teresa 
Lemos, Rogério Pinto 


CONSELHO REDACTORIAL 


A. Ferreira dos Santos, À. Romão 
Dias, Aquiles Gomes, Carlos Ros- 
maninho, Domingos Moura, F. Sou- 
sa Marques, H. Nunes Garcia, J. 
Borges da Silva, J. P. Figanier, L. 
Aires-Barros, L. Borges Teixeira 


DIRECÇÃO, REDACÇÃO 
E ADMINISTRAÇÃO 


Av. Rovisco Pais, 1.S. T. — Lisboa 
Telefone 731549 
PROPRIETÁRIO E EDITOR 
Ma E ts ds 
* 
ASSINATURAS: BA” 100º 
Continente e Ilhas 95800 1808. 
Ultramar, Brasil 
Espanha ...... 100800 71908. 
Estrangeiro ........ -— 2008. 
Número avulso ... = 2558. 


Não se publica em Agosto e 
Setembro 


Os artigos assinados são da 
exclusiva responsabilidade dos au- 
tores. 


* 


COMPOSIÇÃO E IMPRESSÃO NA 


MOVIGRÁFICA 
R. do Espírito Santo, 19-A Odivelas 


DEZEMBRO DE 1971 


| 


ANO XLVI 


VOLUME XXXIV PUBLICAÇÃO MENSAL 


SUMÁRIO 


117 — 


o) — 


143 — 


153 — 


165 — 


TOU = 


tis = 


CAPA — 


FERRY BORGES -— Dynamic loads (in par- 
ticular wind and earthquake loads). 
Solicitações dinâmicas (em particular causadas por 
ventos e sismos). 


DUARTE RAMOS -— Significado do campo 
eléctrico no arvorejamento de polímeros devido 
a uma heterogeneidade elíptica e hiperbólica. 
The importance of electrical field distribuition by 
treenig of polymers with an ellipsoidal hyperbolical 
electrode. 


PERES RODRIGUES — Anisotropia das rochas- 
-superfícies mais prováveis das tensões de 
rotura e dos módulos de elasticidade. 
Anisotropy of Rocks. More probable surfaces of ulti- 
mate stresses and moduli of elasticity. 


M. BARROSO — Caldas de injecção à base 
de cimento. 
Injection mixtures based on cement. 


SALGADO PRATA — Perspectivas da evolu- 
ção dos critérios de aceitação dos defeitos de 
compacidade em soldaduras. 
Evolution perpectives of criteria of acceptance of 
compacity lacks in welder joints. 


VIDA ESCOLAR E ASSOCIATIVA. 
Valadares Tavares e alunos. — Uma expe- 
riência no ensino. 


BIBLIOGRAFIA. 


B. N. R. 1. — 80. Uma das máquinas utilizadas pela CP na 
renovação da via. 

A B. N. R. | — 80 executa simultaneamente o ataque, o 
nivelamento e o alinhamento da via. 


INDICE ALFABÉTICO DOS ANUNCIANTES 


Pág. 
Aguiar & Melo, Lda. . ..... XXVI 
ABRE mansa us us é XXV 
Cimento DOC, «mma nr wess 2º da capa 
Companhia de Cimentos Brancos VI 
Construções Técnicas . . .... XVII 
Construções Schréder . . .... XIX 
DOE as cwmskhs sã e diê XI 
Eficex Klenzle . . ....... IV 
Electro-Arco, Lda. ...... EI 
Empresa de Cimentos de Leiria . X 
Empresa Electro Cerâmica ... XIV 
Empresa de Sondagens e Funda- 
ções Teixeira Duarte, Lda. . 3.º da capa 
Empresa Técnica de Equipamen- 
tos EMCCUACÕE sn css ns: XI 
Empresa Ultramarina de Sonda- 
gens e Fundações. . .... XVI 
MM es emas ca ga XVI 
Fundações Franki, Lda. . .... I 
Guedes & Almeida, Lda. .... IV e VIII 
Harker Sumner . ....... XXIV 
MO esa mama es Es 2.* da capa 
DDD. é is Ev KR EMA & XX 
Johann Keller ........ VI 


Pág. 
J. F. Azevedo e Silva. . . .... VII 
E A emu dn E CER d XX 
SOMAR, mui sema Ea XIV 
LIVIOS LOCHICOS . « es ums su V 
MR emma RN ENS XIII 
Mendes de Almeida . ....... XXI 
Monteiro Gomes . ....c.c.. XXV 
NOVOLTA, ADA. «cu cmmma os 2.º da capa 
Quinzena Técnica ...... II 
REMO. anzcsswvama El XXV 
Salvador Caetano . ....... XXIIJ 
DOM sm ssmu nv Es 
Bem pspmsgasmra 4º da capa 
Sociedade Portuguesa Cavan .. XX 
Ds BM. cms Kva gm a iz XXII 
Sondagens e Fundações A. Cavaco V 
Sondagens Ródio, Lda. ..... VI 
MOMO ;xsz wa cwzgi cê II 
Standard Eléctrica . ..... IX 
Vad. DÃO: : ss se pis XV 


Wild Portugal, Lda. . ..... es: 


espira 


ha base 
lia CONStru 


ae 
Eli 
Riu 
Est 
HE 


tão 
ne 


ão y á ERES E 
Et o 


"a 
E 


E eae > | do ARA 
ope o 
si a ce cansa 


To adega 


A elevada gama de diâmetros das estacas Franki per- 
mite a melhor adaptação a cada problema concreto 
de fundações, desde a construção mais ligeira até à 
de grande concentração de cargas muito elevadas. 
Material permanentemente actualizado e pessoal com 
especialização constantemente aprofundada garantem 


a segurança, rapidez e economia das fundações 
executadas. 


Consulte-nos: FUNDAÇÕES FRANKI, LDA. — PORTUGAL 


TECNICA | 


material de 
trabalhos públicos 

equipamento eléctrico 

climatização 


especialistas e compradores 
este guia é-vos destinado 


A Quinzena Técnica de Paris que agrupa : EXPOMAT (25 de Maio - 4 de Junho de 1972) 
oitavo salão internacional de materiais de Trabalhos Públicos e de construção - 971 exposi- 
tores de 18 nações. BIENAL DO EQUIPAMENTO ELECTRICO E DA AUTOMATIZAÇÃO DOS 
PROCESSOS INDUSTRIAIS (31 de Maio - B de Junho de 1972) 494 expositores de 18 
nações. INTERCLIMA (2 - 8 de Junho de 1972) : aquecimento, retrigeração e climatização - 

450 expositores de 20 nações. ...scaba de publicar para vós um Guia da exposição. Com 
informações pormenorizadas sobre os salões, Inclui um cartão para um programa de visitas 
às obras de “Paris - 1975": Bairro Délense, Adroporto de Roissy em França, Periféricos, torre 
Maina-Montparnasse... pontos de observação previlegiados para os técnicos do mundo Inteiro. 


“Quinzena 
Tica — 


25 DE ” 
= das JUNHO DE 1972 


 dristaeias «sde Paris 

Di acres ss RE E ELO a Do sad C ES bo EPE O NON TEL ONO Viet as a Dona A E ad - 
k t 
4 | 
E EMPRESA q 
o. t 
MORADA - 
4 deseja rátis 0 da 0x una Técnica de Paris ” destinada sos técnicos ' 


uia 
Ê poriqueses mara AMBRE De COMMERCE PRANÇARE DU PORTUGAL. Av. renal "ês Chaves 8 
8 41. 6º Dio LISBOA, 4-97-48 


TFCNICA II 


fundações 
sondagens 


Sopecate 


especialista hã $1/4% de século 


— —e— CS 


rua do arsenal 146-2º tit 34010 +» Lisboa 


ELECTRO-ARCO, L.' 


ELÉCTRODOS + POSTOS + ACESSÓRIOS 
PARA A SOLDADURA ELÉCTRICA 


“* 


MATERIAL APROVADO PELOS 


LLOYD'S REGISTER OF SHIPPING 
BUREAU VERITAS 
GERMANISCHER LLOYD 


LISBOA 
Rua Silva Carvalho, 239 
Tels. 683649-684893 
FÁBRICA E LABORATÓRIOS 
Venda Nova — Amadora 
Tels. 9706 13-971733 
PORTO 
Rua do Bolhão, 216 
Tel. 21277 - 394473 


<gora 
para distâncias 


até 2000m 


WILD 
DISTOMAT 
DI 10 


O mais recente modelo do único aparelho de medição 
de distâncias, por infra-vermelhos, em serviço desde 


3 anos, permite alcance até 2000 m 
e é adaptável aos teodolitos WILD TIA, Ti6 e T2. 


Esta combinação é a ideal para poligonações, 
levantamentos por radiação e para implantações. 


Economiza tempo, trabalho, dinheiro 
Aumenta precisão, rendimento, lucro 
Aprovado desde 3 anos em todá a parte do mundo 


Garantido pelos nossos serviços técnicos 


MEERSBRUOS OG 


Solicitem prospecto detalhado ou demonstração 
aos representantes exclusivos 


WILD PORTUGAL, LDA. 


PRAÇA DAS ÁGUAS LIVRES, 8 * TEL. 689112 x  LISBOA-2 
TECNICA HI 


ESCAVADORAS HIDRÁULICAS 
S/RODAS -—- S/LAGARTAS 


— ENTREGAS IMED 


—— e— 


IATAS 


qual: RCE » 


4 À - o . 
as aê » 
+44 + ot o o 
fo O + 
o Ro +. E cent da 2 = 
E a, a o a. a E a > . 
»* o. m Ea o E a ” . J q nx : gp" ço 


DISTRIBUIDOR EXCLUSIVO 


GUEDES & ALMEIDA, LDA. 


RUA ÁUREA, 181-2.º — TELEFONES 325011"327845 —- LISBOA 


- - instantâneamente, 
silenciosamente, 
com economia e 
segurança. 


-- «programa de 
calculadores electrónicos 
de CASIO COMPUTER, Co. LTD 
TÓKYO JAPÃO 


- « «diversos modelos com 
as mais evoluidas 
características para toda a 
espécie de cálculos. 


Demonstrações sem compromisso 


«12 Modelos 
diferentes. 
- 1 Programável. 


REPRESENTANTES EXCLUSIVOS: 


EFICEX KIENZLE LISBOA: Av. João XXI, 4A-tels 727028-727019 
DISTRIBUIDORES NO NORTE: 


HORBE.LDA] porto: Rua Faria Guimarães, 105- tel 433.46 


TECNICA |V 


Dimensionamento de Secções 
de Betão Armado em Relação 


à Rolura 


Matérias 


— Flexão simples de secções rectangulares 
simplesmente armadas 

— Flexão simples de secções rectangulares 
duplamente armadas 
a) com armadura de compressão reduzida 
b) com armadura de tracção reduzida 

— Flexão simples de secções T duplamente 
armadas 

— Flexão simples desecções T duplamente 
armadas 

— (Compressão simples 

— Flexão composta de secções rectangulares 
a) com armadura simétrica 
b) com armadura de tracção reduzida 
c) com armadura de compressão reduzida 

— Flexão composta desviada de secções rec- 
tangulares 

— Flexão composta de secções em T 
a) com armadura de tracção reduzida 
b) comarmadura de compressão reduzida 

— Flexão simples e flexão composta de sec- 
ções circulares cheias 

— Flexão simples e flexão composta de sec 
ções circulares ocas 

— Esforço transverso 


Preço 150$00 


Dimensionamento de Secções 
Circulares de Betão Armado 


em Relação à Roltura 
Matérias 


— Flexão simples e flexão composta de sec- 
ções circulares cheias 

— Flexão simples e flexão composta de sec- 
ções circulares ocas 


Preço 60$00 


Por J. M. MADEIRA COSTA 
ENG.º CIVIL 


À venda na TÉCNICA e nas LIVRARIAS 


Água: 
ndo 


e Pesquisas e captações de água 
e Estudos hidrogeológicos 

e Sondagens e Fundações 

e Rebaixamento de níveis aquíferos 


ACAVACO 


R Rodrigo da Fonseca, 62, r. c-Esq. - Ltsboa-1 
Tel. 56 1171/2/3/4 


TECNICA MV 


"| O EMPREGO DO CIMENTO BRANCO... 


permite acabamentos 
mais perfeitos, mais 
duradouros e muito 
mais económicos 


-stude a vantagem do emprego do 


CIMENTO BRANCO LUSO 


IBRA 


No Sul: SCIAL — T. do Corpo Santo, 15 — Telef. 32 7377 — Lisboa 
No Norte: SCIAL — R. António Feliciano de Castilho, 1126 — Telef. 97 09 32 — Areosa — Porto 


SONDAGENS RÓDIO, LD. 


R.S. Bento, 644- 3" 
LISBOA 


TELEFONE: 68 80 96 


“a ” SONDAGENS GEOLÓGICAS 


Pressão 


ESTUDOS GEOTÉCNICOS 
CONSOLIDAÇÃO DE FUNDAÇÕES 


instalações de 

Hidráulica Horizontal 
para introdução de Ê PAREDES CONTÍNUAS 
passagens inforiores e canalizações | IMPERMBABILIZAÇÃO DE OBRAS 
sem interrupção do tráfego HIDRÁULICAS 
ferrovióério ou rodoviário Ti CAPTAÇÃO DE AGUAS 


ESTACAS 


JORANNH KELLER 
ELITE 


TECNICA MI 


PROLONGUE A VIDA DAS SUAS 
LÂMPADAS FLUORESCENTES com 


Dstastros PERFEKTSTART ENO 


SEM ARRANCADOR 


Arranque controlado pelo sistema termo-eléctrico 
PERFEKTSTART o qual incorpora um relais térmico 
que proporciona o pré-aquecimento dos cátodos das 
lâmpadas 


po 
E SMA 
+ Va go 
*, 
E 
Sa 


RE 


vinnsas d+ 


nn 


IMPAR 


Balastro PERFEKSTART, em 
corte mostrando o relais 
térmico do dispositivo de 
arranque 


PRINCIPAIS CARACTERÍSTICAS 


— Disponiveis nas potências de 20, 30, 40 e 65W  — Livre escolha do tipo de lâmpada que mais con- 
— Garantidos por*5 anos contra defeitos de fabrico vier, para arranque normal ou rápido 
— Ausência de intermitências ou repetidas tenta- -—- Ausência de indicações especiais para a mon- 


tivas infrutiferas durante o arranque tagem ou de dispositivos auxiliares para melho- 
— Ausência de problemas com arrancadores. Tor- rar o arranque 
na-se desnecessária a sua substituição 


— Lâmpadas fluorescentes de boa qualidade dei- 
— Duração da vida da lâmpada pelo menos dupli- xam de enegrecer nas extremidades 


cada. É possivel contabilizar exactamente a — Ausência de interferências 
economia que se obtém com o seu funciona- rádio ou televisão 


mento e manutenção — Ausência de percas de corrente — bz'xo aque- 
— Funcionamento práticamente silencioso graças cimento 


a novas técnicas de fabrico incorporadas no Marca de qualidade da ASE para utilização em 


núcleo magnético locais secos ou húmidos 


JF. DE AZEVEDO E SILVA & CL: 


Trv. Fábrica das Sedas,22-A + Tel. 6593711-654165 + LISBOA 


nos receptores de 


"TECNICA VII 


WAaBCO 
tg PP MOTO-SORAPERS AUTO-GARREGADORES 


CAPACIDADES ATÉ 34 JARDAS CUBICAS 


NÃO NECESSITA DE TRACTOR PARA EMPURRE! 
UMA SÓ MÁQUINA... UM SÓ OPERADOR... 


DISTRIBUIDOR EXCLUSIVO 


GUEDAL - Guedes & Almeida, Lda. 


R. ÁUREA, 181-2.º º LISBOA-2 ” TEL. 32 50 80 


R. ALFREDO CUNHA, 538 2 MATOSINHOS e TEL. 93 45 62 
TECNICA VII 


À maior rede de especialistas 
“na Europa antecipa-se 


aos seus problemas | 


| O Grupo de Componentes da ITT 
Europa tem uma ambição: 
Não permitir a interrupção 

| da actividade dos seus 

| clientes por falta de 

| componentes. 

| Conhecendo em todo o 

| pormenor as caracteristicas 

dos mercados podemos 

| antecipar-nos à resolução dos 

seus problemas, Por isso temos: 

na Europa Ocidental 4 modernos 

centros de investigação científica | 

que apoiam as divisões comerciais. À 


Somos responsáveis por alguns resultados 
impressionantes: O primeiro 
estabilizador cardiaco baseado na 
tecnologia que a ITT desenvolveu 
para os seus circuitos 
integrados; a produção 
de amplificadores 
para cabos 
submarinos 

que, submersos 
a grandes 


profundidades, têm 

20 anos; reduzimos para 3 

Natureza cerca de 3 milhões de anos 
fenômenos fisicos, base do 


Psi o) 
que funcionar com 
semanas o ciclo de produção de 
a conseguir. 
Os nossos programas de pesquisa 
baseiam-se essencialmente numa 
continua procura de novas 
desenvolvimento de 
If | 


absoluta segurança durante 
cristais para rádio o que levaria à 
matérias-primas e de novos 
técnicas e processos que 


| levam ao aparecimento 

Exemplos de componentes ITT | co Hovor compoNarada 
a Sao 

São estes motores e ventiladores largamente | 

usados em aplicações tais como: máquinas pisca aaa 

de executivo, gravadores de fita magnética, 

aparelhagem doméstica, projectores, 

aquecedores de água, sistemas de ventilação 

e refrigeração, etc 

A gama de motores ITT está prevista para 

operação unifásica ou trifásica com potências 

nominais de 0,7 a 78 W. As potências dos 

motores de corrente continua vão desde 0,25 

a 25 W, Existem ainda modelos acoplados 

a redutores de velocidade. 


Para informações completas dirija-se a: 
STANDARD ELÉCTRICA, S.A. R.L. 

Divisão ITT de Vendas de Componentes 

Av. da República, 60, 2.º-D. 

Lisboa 1 

Telef. 765491 | COMPONENTES 


TECNICA IX 


E E 
| 


EMPRESA DE CIMENTOS 
DE LEIRIA 


ED» Bio ER o ddô 


CIMENTO SUPER-VULCANO 
CALCINA 


CAL HIDRAULICA 


SÍMBOLO DE QUALIDADE 


Sede: Rua Braamcamp, 7 LISBOA ; Fábrica: Maceira — Liz 


TECNICA X 


Comunicações 


equipamentos 
e sistemas completos 


+ Cuidadosamente 
projectados 


- Amplamente 
experimentados 


* Completamente 
integrados 


THE TECHNICAL MATERIEL CORPORATION 


Registador Sanborm 


HEWLETT PACKARD 


Registador digital 


aparelhagem electrónica Dymec 
de medida e registo 
Osciioscópios, Oscilógrafos, Voltímetros 
electrónicos, Geradores de sinal, 


Contadores electrónicos, Registadores 
digitais, Fontes de alimentação 


Oscilocónio HP 


Contador EHectrónico HP 


VINTA:, A EMPRESA TECNICA DE EQUIPAMENTOS FLECTRICOS, SARL 


TECNICA XI 


Motores eléctricos 


E AE 


a». 


S 


* para todos 
os fins 
* industriais 


E a Potências até 3000 CV 


FFAGEC, a maior organização do país no ramo electrotécnico 


ras 


e 


TECNICA XI 


ANO XLVI 


| 


VOLUME XXXIV PUBLICAÇÃO MENSAL 


C. D. U. 624.042.8 


DYNAMIC LOADS (IN PARTICULAR WIND AND 
EARTHQUAKE LOADS) * 


NUMERO 413 DEZEMBRO DE 1971 
SUMMARY 
Present information on dynamic loads, such as 


those due to wind, earthquake, traffic, machinery and 
blast, in briefly reviewed. 

Attention is focused on the first two of the men- 
tioned loads. The main problems considered are: i) mea- 
surement of acting loads, and ii) their analytical and ex- 
perimental idealization. 

The general lines along which research is to be 
continued are suggested as a basis for discussion. 


1 — INTRODUCTION 


Progress in the structural field implies the impro- 
vement of knowledge of both the acting loads and the 
structural behaviour under these loads. In general it can 
bc said that the actual knowledge of dynamic loads is 
not great. So, it is easy to understand the interest of not 
only obtaining data that allow a better definition of these 
loads, but also of establishing the most convenient 
methodology for their analysis. 

For the distinction between static and dynamic 
loads a first difficulty arises from the fact that this dis- 
tinction cannot be established on the basis of the loads 
themselves but has to consider the type of structure on 
which the loads act. In fact, a load has to be conside- 
red as dynamic when its variability in time is such that 
to study the structural behaviour the effects of the iner- 
tia fo the masses of the structure cannot be neglected. 
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RESUMO 


Acrescenta-se uma sintese da informação actual- 
mente disponivel relativa às principais solicitações dinã- 
micas que actuam em estruturas de engenharia civil, 
tais como vento, sismos, tráfego, vibração de máquinas 
e explosões. 

Presta-se particular atenção às solicitações dos dois 
primeiros tipos mencionados e, para estas, consideram- 
-se os problemas principais de: i) medição das solicita- 
ções actuantes: ii) sua idealização analítica e experi- 
mental. 

Como conclusão e tema de discussão, propôem-se 
linhas mestras para prosseguimento da investigação 
neste domínio. 


A second difficulty concerning the definition of 
dynamic loads arises from the interaction between the 
structural behaviour and the loads themselves. Conside- 
ring the dynamic loads as a system of forces variable 
in time, these forces are often directly influenced by 
the behaviour of the structure. This for instance occurs 
in wind actions due to aeroelastic effects and in earth 
quake problems owing to the interaction between the 
structure and the soil. In the present report it is tried to 
define the loads in a fundamental way, in order that 
these interactions may be studied by considering the 
total behaviour of the overall system. 

To limit the scope of the report attention is focus- 
sed on the loads that act on civil engineering structures 
such as buildings, bridges and towers. 

The types of loads to be referred are those due to: 
wind, earthquake, trafíic, machinery and blast. The two 


* Introductory report presented at the VIll Congress of the International Association for Bridge and Structural Engineering, New 
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first types of loads which in general are more important 
are studied in more detail, 

For each type of load the main problems to be con- 
sidered are: i) measurement of the acting loads, ii) their 
analytical and experimental representations. The dis- 
cussion of the most convenient methods for cealing with 
these problems is very important and, as referred, it is 
hoped that during the Congress the discussion shall co- 
ver methodology and not be limited to the presentation 
of results. 

For the purpose of choice of method two funda- 
mental types of representation are considered: deter- 
miniístic and stochastic. The deterministic representation 
implies the complete knowledge of the variation of the 
loads with time. The stochastic representation only im- 
plies the knowledge of the statistical distributions of 
these variations. 

Machine vibrations or blast pressures may be exam- 
ples of loads for which a deterministic approach applies. 
Wind and earthquake loads in general can only be con- 
veniently represented by stochastic schemes. 


The improvement of load knowledge implies that 
systematic measurements are performed in nature. The 
dynamics character of the loads makes these measu- 
rements particularly difficult. The difficulties are of two 
principal kinds. The first comes from the fact that often 
it is not possible or practical to measure directly the 
acting forces but only other magnitudes from which the 
forces are to be derived. That is for instance the case for 
the wind for which it is the velocity and not the pressure 
that has to be considered. So a new problem has to be 
solved: to transform the velocities in pressures. 


The second difficulty comes from the response of 
the measuring device in function of the frequencies in- 
volved. In fact the measuring devices behave in general 
as low-pass filters, have often a non-flat frequency re- 
sponse and even in some cases are non-linear, that is, 
the response is a function of the mean intensity of the 
load. For a correct interpretation of the results, the 
measuring devices have to be considered as transfer sys- 
tems and the characteristics of their response accurate- 
ly studied. 


Difficulties of this type appear for instance in rela- 
tion to the use of anemometers. The results available 
have to be interpreted by duly taking into account the 
dynamic characteristics of the apparatus. 


Once the nature of the loads is well understood 
the convenient representations of the loads to be used 
in analytical or experimental studies have to be dis- 
cussed. The representations must be of a nature as fun- 
damental as possible in order that they can be applied 
with generality. Also the aralytical representations of the 
Icods must fit the general theories to be used for study- 
ing structural behaviour. Particularly in the stochastic 
scheme the fundamental concepts of the theory of random 
vibrations [1,2] must be respected. 

It is also convenient that analytical representations 
are as simple as possible. The influence of the intro- 
duced simplifications must be analysed and the range 
in which they apply must be defined. 
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In the case of permanent loads the definition of 
duration is not important. In the case of non-perma- 
nent loads it is convenient to split up them by the 
consideration of intervals with given durations. So wind 
and earthquake loads shall be represented by time series 
having each a given duration. A further usual simpli- 
fication consists in supposing that the statistical distri- 
bution from which the time series derives does not 
change along the time for the assumed duration; tha! 
is to suppose that the phenomena are stationary. 

The experimental representation of the dynamic 
loads must also satisfy some general conditions. Loads 
representing the actual ones are to be imposed on mod- 
eis and behaviour of these models has to be inter- 
preted using theories of similitude. The frame in which 
similitude applies imposes conditions on the lvad repre- 
sentation. When interpreting the results the limits within 
which the tests were performed have to be duly con- 
sidered. 

So, for instance, wind tunnel tests in general dis- 
regard wind gradients and do not respect the similitude 
for turbulence. Dynamic tests for the study of earthquake 
problems often use vibrations that are far from represen- 
ting earthquake movements. 

To allow a comparison between analytical and 
experimental results it is convenient that analytical and 
experimental representations of the loads are in accor- 
dance. 

Finally, the problem of load forecastiny has to be 
considered. This forecasting must in general be estab- 
lished on statistical bases. Statistical concepts of safety 
may then be applied. 

Is must be well understood that this statistical 
forecasting has nothing to do with the (deterministic or 
stochastic) scheme adopted for the representation of the 
load variations in time. 

The general problems concerning dynamic loads on 
structures have been recently studied in several symposia, 
among which the «RILEM Symposium on the Measure- 
ment of Dynamic Effects and Vibrations of Construc- 
tions» held in Budapest in 1963, the «Symposium on 
Vibration in Civil Engineering» held in London in 1965 
[3], and the «RILEM Symposium on the Effects of 
Repeated Loading of Materials and Structures» held ir 
Mexico in 1966. The particular problems concerning 
wind and earthquake loads have been dealt with in 
special meetings that are referred below. Davenport in 
the report presented at the symposium held in Mexico 
[4] makes an interesting general analysis and a com- 
parison between wind and earthquake loads. 


2 — WIND LOADS 

2.1 — Nature of wind 

A good understanding of wind phenomena implies 
the general knowledge of wind causes. 

As is well known, the wind velocity increases 


with the height above ground and for heights of about 
300 to 600 m reaches a limit value (gradient velocity) 
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that is mainly governed by pressure gradients (direct 
connected with thermal effects) and geostrophic accel- 
erations. This gradient velocity can be analytically related 
to the mentioned causes. 

The wind velocity rapidly varies with time and from 
point to point. This variation being rancom the flow is 
turbulent. 

For the study of wind loads it is convenient to con- 
sider intervals of time with a given duration (for instance 
10 minutes or 1 hour) and to compute for these time 
intervals the correspondent mean wind velocity. The 
rapidly varying velocity component then corresponds to 
turbulence. 


Recent studies on the structure of wind that duly 
consider turbulence have been principally performed in 
connection with air pollution [5] and aeronautics [6]. 
Much information coming from these sources is now 
available but only a small part can be used for studying 
the wind action on structures. Even so, it must be re- 
cognized that the methodology used in these studies is 
the convenient one. Modern studies on the wind action 
on structures follow the same lines [7, 8]. 


2.2 — Wind measurement 


The measurement of the wind velocity may be divi- 
ded into two different problems. The first concerns the 
measurement of the low frequency component that 
corresponds with the mean velocity. The second con 
cerns the high frequency component due to turbulence. 


Practically all over the world the measurement of 
mean wind velocities has been in charge of the meteo- 
rological services and a large information is available on 
them. On the contrary measurements of turbulence have 
been performed in relation with particular researches 
only. 
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The meteorological services also indicate maximum 
velocities but these are difficult to interpret and for the 
time being are unreliable for use in structural design. 

The different types of anemometers for measuring 
wind velocity can be classified according to their range 
of frequency. 

Among the low frequency types are: Pitot static 
tubes, propeller and vane anemometers. These instru- 
ments have in general cut-off frequencies of about 
0.1 c/s. 

In the medium range of frequency are membrane 
or vibrating mass anemometers with cut-off frequencies 
of the order of 20 c/s. 

Finally in the high range of frequency can be 
considered the hot wire and electric discharge anemo- 
meters that respond up to frequencies of about 1000 c'/s. 


Fig. 1 presents a typical record obtained using a 
Pitot tube anemometer. 


For the interpretation of the records of these dif- 
ferent types of anemometers it is necessary not only to 
calibrate them statically in wind tunnels but also to 
determino their transfer functions. Besides, their behav- 
iour being ofíten not perfectly linear, the transfer func- 
tions have to be determined for different values of the 
static component. 


Measurements in turbulent flows, and particularly 
the measurement of turbulence itself, involve important 
difficulties [9]. These difficulties are much reduced if 
only mean wind velocities are considered. Owing to the 
presented reasons, only the mean values of wind velocity 
constitute a set of valuable data from which the maxi- 
mum velocities can be derived, The relashionship bet- 
ween mean and maximum velocities has to be based on 
the actual knowledge of the structure of wind. The 
spreading of the correct measurement in the high fre- 
quency range is most desirable. The data then obtained 
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shall allow not only a better understanding of turbulence 
but also direct estimates of maximum velocities. 


2.3 — Representation of wind 
2.3.1 — Mean wind velocity 


For the study of the wind loads on a structure it 
only interests to consider a limited spatial domain sur- 
rounding the structure. Also, as referred, only a given 
interval of time has to be considered, for instance 10 
minutes or 1 hour. 

In general it shall be reasonable to suppose that in 
tho given domain the mean velocity is represented by a 
horizontal component that only varies in function of the 
height. To define the field of mean wind velocities it is 
then necessary to define the law of variation of the mean 
wind velocity with height only. 

For the range of height that interests civil engi- 
neering and for winds of high velocity (those of interest 
for the design of structures) the mean velocity at level z, 


U, , may be expressed by a law of the type: 


cais : NU. 
tj a 


where Uc is the gradient velocity reached at level 26: 
The values of Zs and 1/z cepend on the roughness 


of the ground, 


Às the wind velocities are in general measured at 
heights of about 10 m above the ground it is convenient 


to use for reference the velocity at this height, Uis 
and not the gradient velocity. 

Fig. 2 gives the mean wind velocity profiles, acor- 
ding to Davenport [8] for three typical conditions of 
ground roughness. 
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Fig. 2 — Mean wind velocity profiles according to Davenport (8) 
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To each location there corresponds a statistical dis- 


tribution of Qi , that defines the probability of this ve- 
locity being exceeded during a given interval of time. 


2.3.2 — Turbulence of wind 


Due to the existence of turbulence it is necessary 
to add to the mean velocity a varying velocity defined 
by the components u, v, w (respectively longitudinal, 
transverse and vertical). 

These components are supposed to have a random 
variation in time and space. The statistical definition of 
these variations can be performed only by assuming 
simplifying hypotheses. 

A first reasonable hypothesis consists in supposing 
that the turbulence is stationary during the time interval 
considered. The variation of velocity at a point may then 
be statistically described by the variation of the spectral 
densities of the velocities in function of the frequency 
or by the Fourier transform of these spectral densities, 
the autocorrelations. Assuming homogeneity at the dif- 
ferent levels, the spectral densities of velocity may then 
vary in function of the mean wind velocity and the 
height above ground. 

Spectral densities in time are not sufficient to define 
turbulence completely. It is also necessary to consider 
tho statistical variation in space cescribed by the space 
spectral densities or space correlations. 

For aeronautical problems the turbulence of the ver- 
tical component is very important [6]. The same does 
not apply, in general, to the wind actions on structures. 
In this case the longitudinal component is of particular 
interest. 

Fig. 3 presents the spectral density of the longitu- 
dinal component of wind velocity proposed by Daven- 
port [8]. It represents the mean value of measurements 
performed at heights ranging from 8 to 150 m. As means 
were taken, this implies to consider that turbulence does 
not vary with height in the range considered. 


For representing the spectra obtained in different 
concitions by a single curve, Davenport divides the 
spectral censity S(f) by a coeficient K, a surface drag 
coeificient depending on the ground roughness, and by 


TEA the square value of the mean velocity at the refe- 
rence level. The abscissae represent in a logarithmic 
scale not cirectly the frequencies but wave numbers or 
wave lenghts. The wave mumber, expressed in cycles 
per meter, may be reduced to the frequency (c/s) by 
multiplying it by the mean velocity (m/s). The wave 
length is obtained by dividing the mean velocity by the 
frequency. 


Taylor hypothesis consists in supposing the validity 


of the transformation x = Ut, that is, it establishes an 
equivalence between the variations in space (x) and 
in time (t). According to this hypothesis the presented 
spectral density in time may also considered as a spec- 
tral density in the longitudinal direction. 

As the abscissae are represented in a logarithmic 
scale it is convenient to multiply the ordinates by the 
frequency f in order that the integral of the spectral den- 
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Fig. 3 — Spectral density of longitudinal component of wind velocity according to Davemport (8) 


sity represents the mean square value of the velocity 
fluctuations. 

For wave lengths smaller than 500 m the proposed 
spectrum fits Kolmogorov's law [10]. This law relates 
the spectral density S (2) to the wave number, >, by 
an expression of the type S (2) = 54 A 

To get quantitative information that may be used 
in structural design the spectrum of fig. 3 is plotted in 
fig. 4 for a mean velocity of wind of 20 m/s. 

This figure shows that the reduced spectral density, 
f S(Z) . 
——— |, has a maximum between 1 and 2 cycles/mi- 
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Fig. 4 — Reduced spectral densities of wind velocity and earth- 
quake acceleration 
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nute. As the frequency increases the spectral density ra- 
pidly decreases and for practical purposes may be con- 
sidered to vanish above 1 or 2 c/s. 


The relative intensity of turbulence that corresponds 


] 
to the spectra of fig. 3 and 4 is given by = 
10 
For K = 0.005 this corresponds to a value of 0.17 
that is in accordance with the usual intensities of tur- 
bulence, between 10 and 20%, indicated by Pasquill 
[10]. 

The variation of the longitudinal component of 
wind velocity in transverse directions can only be de- 
fined considering the correlation in these directions or 
the cross spectral densities S (AI, f) for points at dif- 
ferent distances, A 


=V6K. 


Davenport [11] relates the cross spectral densities 
to the spectral density in time, S (f), by means of a 
coefficient 


S (al, +) 
R(441) = E = 
0 Ca 


where c is a constant. = has the dimensions of a 
c 


length and can be interpreted as a correlation scale. 


The distance at which wind velocities are correlated 
is therefore inversely proportional to the frequency and 
to the coefficient c. Table | indicates for stable atmos- 
pheric conditions (that are those of interest in the pre- 
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sent case) the order of magnitude of the value of c for 
different directions and turbulence components [11]. 


Table | 
Direction | Turbulence é 
of AÍ component 
longitudinal 8 
longitudinal Ea = 
transverse 6 
longitudinal 40 
transverse di ! 
transverse 25 
longitudinal E) 
vertical RSRS 
transverse T 


The fact of the correlation scales being smaller in 
the transverse than in the longitudinal direction indicates 
that for strong winds eddies are elongated. 


2.4 — Wind actions 


The knowledge of wind velocities is not sufficient 
to define wind actions. It is necessary to know how to 
transform the velocities in pressures. Duly considering 
the dynamic character of the phenomena this problem is 
not yet satisfactorily solved. 

Supposing the wind to be a uniform flow, much in- 
formation is available that allows to transform the wind 
velocity into local or total pressures. This information 
is collected in wind codes under the form of pressure 
coefficients. 

Jensen [12] has shown that if the uniform flow is 
substituted by a boundary layer profile, as indicated in 
fig. 2, important variations of the pressure coefficients 
are obtained. 

Considering the turbulent character of wind the 
problem is much more involved. In fact to the turbulence 
inherent to the wind it is necessary to add the turbulence 
created by the structure. 

The turbulence created by the structure derives 
mainly from vortex excitation and may be independent 
or dependent of the deformability of the structure itself. 

The first type of phenomena corresponds with von 
Kármán eddies. The frequency, f, of the vortex-shedding 


is related to the mean wind velocity, U, by the Strouhal 


fD 
number, S ME ci , where D is a typical length. As 


there is a dominant frequency f this is not a true tur- 
bulence; it is often called a quasi-turbulence. A review 
of the problem of vortex-shedding for rigid circular cylin- 
ders is presented by Lienhard [13]. 

The excitation of vibration due to the deformability 
of the structures is a still more complex phenomenon. A 
simple case where the mechanism of the vibrations can 
be easily understood is the following. 
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Consider a horizontal wind and a structure vibrating 
vertically. The vibration of the structure corresponds 
to a transverse component of velocity. To combine 
this transverse component with the longitudinal one is 
equivalent to consider an oblique incidence. If the pro- 
file of the body is such that for this oblique incidence 
there is a negative lift, this lift force tends to increase 
the vibration of the system and may be considered a 
negative damping. If this negative damping exceeds the 
positive damping of the structure, self-excited vibrations 
of increasing magnitude will occur. 


For systems with several degrees of freedom the 
phenomena of aerodynamic instability may be associated 
with coupling of different modes giving rise to flutter 
phenomena. 


Finallv it may also occur that the turbulence pro- 
duced by one structure influences neighbouring ones. 


Scruton [14] reviews the different aspects of 
wind-excited vibrations of structures. Owing to the com- 
plexity of the involved phenomena it can be said that 
general analytical solutions cannot yet be obtained. 
Model tests in wind tunnels yield very useful result. 


2.5 — Model tests 


Usually, when building a wind tunnel it is sought 
to reduce turbulence as much as possible in order to 
obtain a uniform air flow. A tunnel in these conditions is 
well adapted for determining pressure coefficients, but 
obviously does not allow to study the influence of the 
turbulence of natural wind. 

Recently several tunnelis were built [8] to repro- 
duce the variation of wind velocity with height. This is 
obtained by increasing the roughness of the deck wall in 
comparison with that of the other walls. Fig. 5 repre- 
sents the boundary layer wind tunnel of the University 
of Western Ontario. 

In future it would be desirable to be able accurately 
to control the turbulence of wind tunnels in order to 
perform dynamic tests with the turbulence to scale. 


2.6 — Data for design 


Among the data for design it is necessary to deiine 
ths maximum wind velocities to which the pressure 
coefficients apply. The maximum velocities at present 
used are not currently defined on a statistical basis. 

lt is important to be able to relate maximum and 
mean velocities. Several authors studied this problem 
[15, 16]. Castanheta [17] adopting the velocity power 
spectrum indicated by Davenport [8] computed the 
statistic distribution of the ratio of the maximum velo- 
cities, U,. to the mean velocities recorded in time in- 
tervals to = 10 minutes or 1 hour. The maximum ve- 
locities, VU. are those that would be recorded by ideal 
filters with sharp cut-off frequencies f = 0.2, 0.5, 1.0 
and 2,0 c/s. 

Table Il indicates the obtained ratios of the mean 
maximum velocities to the mean velocities. These values 
refer to a height of 10 m above ground. 
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Fig. 5 — Boundary-layer wind tunnel of the University of Western Ontario 


Table Il 
to = 10 minutes | to = 1 hour 
Cut-off ———. 
fre- Roughness Roughness 
FarrE RA, coefficient, K (*) coefficient, K (*) 
fi(c/s)-— — — 


K=0.005 K= =0.010 = SOIS K= K=0.005 k= 0.010 K=0.015 


os [177 1.54 | 1.75 | 1.92 


0.5 | 1.50 | 1.71 | 1.87 | 1.58 | 1.83 | 2.02 


1.0 153 | 175. 1.92 |! 1.63 | 1.88 | 2.07 
2.0 156 | 179) 197 166 | 192 | 213 


(*) As defined by Davenport dl K = 0.005 corresponds 
to a flat ground and 0.015 to rough conditions, as 
those occurring in the centre of a town. 


It is interesting to notice that the variation of the 
coefficients due to the variation of the cut-off frequen- 
cies is practically independent of the roughness coeffi- 
cient and the time interval. 

Taking for reference the frequency of 2.0 c/s reduc- 
tions of 2, 4 and 9% are obtained for the cut-off fre- 
quencies of 1.0, 0.5 and 0.2 c/s, respectively. 

Assuming that a structure behaves like a filter with 
a given cut-off frequency the indicated reductions could 
be applied to the maximum velocities used for design. 
As the pressures are proportional to the square of the 
velocities, the reductions of the pressure would be twice 
the indicated ones. 

In the case of structures with long fronts a further 
reduction can be considered related to the transverse 
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correlation of wind velocities. So far these reductions 
have been mainly established on experimental bases 
[18, 19]. The information now available on the turbu- 
tence of wind is already sufficient lo allow an analytical 
determination of these reductions [17, 20]. 

The above considerations show that reductions of 
wind loads due to the turbulence of wind are only jus- 
tified for structures of very low frequency, such as long 
span suspension bridges, electricity transmission lines 
and high-rise buildings, fig. 4. In fact structures of these 
types present natural frequencies of about 0.1 c/s or 
even less. The same does not apply to ordinary buildings 
(of no more than 10 stories) that have in general fre- 
quencies above 1 c/s. For structures of this type wind 
can be considered as a static load. 

On the other hand for very deformable structures, 
such as suspension bridges, aeroelastic phenomena may 
be paramount importance [211]. 

For improving the knowledge concerning wind loads 
it seems very important to observe the real behaviour of 
structures [22]. Such studies will confirm the assumed 
hypotheses and indicate the most promising research 
lines. 


3 — EARTHQUAKE LOADS 

3.1 — Nature of earthquakes 

It is generally accepted that earthquakes are pro- 
duced by local ruptures of the earth's crust. Earthquake 
vibrations may be felt at very large distances from the 


epicentral zone but, in general, they only affect con- 
structions in a much smaller radius. Within the area in 
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which earthquakes are destructive, soil movements have 
an irregular character as can be appreciated from the 
available records, fig. 6. 

Man-made explosions may also be the cause of soil 
vibrations that in some aspects may be compared with 
the movements due to earthquakes. 

Owing to the earthquake vibrations, water masses 
enter in movement and may produce large waves in the 
sea (tsunami) and also hydrodynamic pressures on struc- 
tures that are submerged in or in contact with water. 
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the number of strong-motion accelerographs is rapidly 
increasing. Even so the most important seismic regions 
are not yet conveniently covered. 


3.3 — Representation of earthquake 


A very important contribution to the representation 
of soil movements due to earthquakes is due to Housner 
[27] who assimilated these movements to a set of ran- 
dom pulses and represented them by acceleration, ve- 
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Fig. 6 — Typical record of earthquake acceleration 


In recent years much effort has been made to in- 
vestigate earthquake causes, to define the seismicity of 
the different regions of the globe, to study the most ade- 
quate methods to measure and represent seismic move- 
ments and to study the behaviour of constructions under 
earthquake loads. Among the abundant literature on these 
subjects, the proceedings of the World Conferences of 
the International Association on Earthquake Engineering 
are of paramount importance [23, 24, 25]. The recent 
effort of UNESCO in coordinating and promoting re- 
search in this field must also be emphasized, 


3.2 — Earthquake measurement 


Although the small amplitude movements of soil 
due to distant earthquakes have been recorded for many 
years, only recently the convenient equipment for recor- 
ding strong motions was developed [26]. 

The apparatus now currently used records the three 
components of the soil acceleration and automatically 
starts the recording when the vertical or one of the hori- 
rizontal components exceed about 0.0] q. 

The dynamic characteristics of this equipment and 
the velocity at which the record is performed are ade- 
quate for further analysis of the seismic vibrations. Also 
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locity or displacement spectra. As defined by Housner, 
the spectrum of a given magnitude, for instance velocity, 
indicates the maximum values of velocity that simple 
oscillators with different natural frequencies and diffe- 
rent damping suffer when subjected to the considered 
motion. Housner computed the spectra of several strong- 
motion accelerograms and proposed to represent earth- 
quake motions by the spectra indicated in fig. 7. 
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Fig. 7 — Spectra of maximum velocities according to Hous- 
ner (27) 
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Another way to describe the randomness of soil 
accelerations consists in defining the power spectral den- 
sity of acceleration. This description uses the fundamen- 
tal magnitudes of the theory of random vibrations that 
were also used to describe wind turbulence. 

The representation adopted by Housner and others 
[28, 29] corresponds to a white noise vibration (vibra- 
tion of constant power spectral density). The variation 
of the power spectral density in function of the fre- 
quency was proposed by Tajimi [30], on basis of preli- 
minary work of Kanai [31], and by Barstein [32]. In a 
study dealing with the probabilistic approach to earth- 
quake-resistant design Rosenblueth [33] reviews the 
different idealizations of seismic loads. 

Bycroft [34] has shown that the velocity spectra 
indicated by Housner were equivalent to a white noise 
vibration limited to the range from 0.2 to 5 c/s, each 
sample having a 30 s duration. He indicates the spectral 
density of 695 cm? s4/c/s as equivalent to the NS 
record of EL Centro 1940 earthquake. 

Recently Ravara [35], Jennings [36] and Arias 
and Petit Laurent [37] used digital computers to deter 
mine direct from the available earthquake records the 
correspondent spectral densities. 

Although the spectra present important fluctuations 
it is possible to deduce a mean law of variation of the 
spectral densities in function of the frequency (fig. 8). 
Analytical expressions for this law have been proposed 
by several authors [30, 36, 38, 39]. It is to be expected 
that the variation of the spectral density will depend on 
the geometry and mechanical properties of the soil. 
Unhappily the information now available does not yet 
permit to quantify this influence nor the one deriving 
from the epicentral distance. 


SPECTRAL DENSITY OF ACCELERATION (cm?s/c/s) 
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Fig. 8 — Variation of the spectral density of acceleration in 


function of frequency 


It must be emphasized that the difference between 
the representation by a white noise of limited range of 
frequencies and the spectrum function corresponding 
to a linear oscillator is not as important as could be 
imagined. In fact mechanical systems always behave as 
filters that cut-off the frequencies above a given limit, 
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and so it is just the same whether loading has a zero 
or non zero spectral density above this limit. Pereira 
[39] compared the response of linear oscillators for 3 
types of spectral functions (fig. 8) and obtained the 
results indicated in fig. 9 that well confirm the above 
conclusion. 
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To represent the soil movements completely it 
would be necessary to define at each point not only 
the variation of acceleration with time but also the cross- 
correlations in different directions. There are at present 
no elements allowing to obtain these data. Even so the 
correlations in time computed by Barstein [32] and 
Arias and Petit Laurent [37] give useful information 
concerning the dimensions of the areas in which the soil 
vibrations can be considered as approximately uniform. 
In fact the time interval for which the correlatio.; falls 
down is of the order of 0.1 s. So even considering high 
velocities of propagation of the seismic waves, points 
at a distance of about 100 m shall have simultaneously 
similar movements. The same may no longer be true 
if the points are more than 500 m apart. 

À representation of earthquakes that allows to study 
the non-uniformity of the vibrations at different points 
was proposed by Bogdanoff, Goldberg and Schiff [40]. 
These authors consider packets of damped oscillatory 
waves with random amplitudes, frequencies, arrival 
times, phases and velocities of propagation. This repre- 
sentation is used to study the longitudinal vibration of a 
suspension bridge. 


3.4 — Structural behaviour 


Analytical methods now available allow to study 
the behaviour of structures under very general hypo- 
theses. Lumped mass methods expressed in matrix form, 
and making use of digital computers, are a powerful tool 
for the dynamic analysis of structures. 


12% 


The earthquake idealization presented above con- 
stitute the necessary basis for performing this dynamic 
analysis. 

When considering earthquake representation it can- 
not be forgotten that, due to the interaction between 
the soil and the foundation, the vibrations in fact applied 
to the structure may be very different form those con- 
sidered. This problem can only be dealt with by duly 
studying the behaviour of the whole structure-soil system 
and assuming that suíficiently for from the structure 
the soil vibration is in accordance with the given repre- 
sentations. 


3.5 — Model tests 


Although the power of analytical methods is nowa- 
days very much increased by the use of digital compu- 
ters many problems cannot yet be conveniently solved 
in this way and for them model studies constitute the 
convenient approach. Both mechanical models and elec- 
tric analogies may be considered as models, although 
electric analogies are in general more close to the ana- 
lytical representation. 

Analog computers have been used with much sucess 
to study dynamic problems concerning earthquake actions 
[14]. Random vibrations may be conveniently studied 
in this way. In fact, random noise generators constitute 


sinusoidal vibrations whose frequency depends on the 
system of springs attached to the table. 

At the Laboratório Nacional de Engenharia Civil, 
model tests have been mainly performed using random 
vibrations, since 1960 [43, 44]. Fig. 10 shows the test 
set-up. The loads applied to the model represent to a 
convenient scale the soil vibration and both the spectrum 
(usually considered a limited range white noise) and the 
duration are reproduced. By performing several tests it is 
possible to Cetermine the mean maximum values of the 
response. By successively increasing the power spectral 
censity it is possible to study the behaviour in the non- 
linear range and even up to rupture. Model studies in 
which random noise vibrations are used have also been 
recently performed in Japan [45, 46]. 

The shakers used in these studies are of the elec- 
tromagnretic type and so only adequate for frequencies 
above 20 c/s. Automatically controlled hydraulic jacx 
systems now produced in some countries are well ade- 
quate to perform model dynamic tests. With these sys- 
tems forces of the order of magnitude of hundreds of 
tons may be applied and these forces can be varied 
up to frequencies of about 20 c/s [47]. As in some 
cases forces can be varied according to any given pro- 
gram, random vibrations may also be applied. It is to be 
expected that systems of this type shall be very useful 
for seismic studies on models. 


Fig. 10 — Random-vibration set-up for model testing 


a standard equipment that may be used to feed the 
analog computer. 

In dynamic tests different techniques can be fol- 
lowed to reproduce soil vibrations [42]. For many years 
sinusoidal shaking tables have been in use with control 
of frequency and amplitude. In other cases vibrations 
have been induced by impacts that produce damped 
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3.6 — Data for design 


In general the codes on earthquake resistant con- 
struction [48] specify seismic coefficients to be applied 
to the weight of the masses existing at the different 
levels, thus allowing to compute the horizontal forces 
for which the structure has to be designed. In this way 
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earthquake dynamic forces are transformed in equivalent 
static ones. 

It is very difficult to establish the seismic coeffi- 
cients in order to contain all the information correspon- 
ding to a complet dynamic analysis. Although rough, 
the design according to seismic coefficients guarantees a 
resistance to horizontal forces, which is very important 
from a practical point of view. 

Considering the difficulties of dynamic analysis, the 
specification of seismic coefficients may be considered 
a convenient procedure to give information to designers, 
mainly for usual types of structures as, for instance, 
ordinary buildings. For structures such as dams, high- 
rise buildings, big bridges and towers it is in general 
necessary to carry out a complete dynamic analysis [49]. 


4 — OTHER TYPES OF DYNAMIC LOADS 


The space available for this preliminairy report only 
allows a very brief reference to the principal types of 
dynamic loads. 

The following considerations are only intended to 
introduce and stimulate the discussion on dynamic loads 
besides those due to wind and earthquakes. 


4.1 — Traffic loads 


The problem of traffic dynamic loads on bridges is 
entirely different according as road or railway bridges are 
concerned. 

For road bridges the actual tendency of codes is to 
distinguish between loads due to traffic congestions and 
loads due to exceptionally heavy vehicles [50]. Interme- 
diate hypotheses do not correspond to extreme condi- 
tions. As traffic congestions may be considered as static 
loads, dynamic behaviour has to be studied only for 
exceptionally heavy vehicles. 

For railway bridges the problem is completely 
different. The repetition of the loads, associated very 
often with inversion of stresses, may produce fatigue. 
The study of dynamic behaviour is of paramount impor- 
tance in this case. 

Another type of dynamic traffic load is the one 
due to landing aircrafts. 

For the improvement of the actual knowledge of 
traffic loads it seems convenient to decompose the dy- 
namic effects into deterministic and stochastic parts. 
These parts would be separately analysed according to 
the respective theories. Thus design rules more accurate 
than the present ones would be obtained. 


4.2 — Machinery loads 


Dynamic loads due to machinery may in general be 
established on a deterministic basis by being assimilated 
to periodic vibrations. In some cases, special types of 
machinery such as ball-mills also produce random vi- 
brations. 

Each machine vibration problem has its own pecu- 
liarities and is is difficult to give general information of 
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interest. Specialized books [51, 52] contain useful in- 
formation on this subject. 


4.3 — Blast loads 


The general term of blast refers both to soil vibra- 
tions and to fluctuations of air pressure due to man-made 
explosions. 

Concerning soil vibrations, blast effects may be 
assimilated to seismic movements. Available records 
show that the accelerations may also be considered as 
random, but the time duration of the vibrations is much 
shorter than the one assumed for earthquakes [53]. 

Studies dealing with quarry blasting [54] relate the 
damage in the constructions and the acceleration or ve- 
locity peaks with the distance to the shot and explosive 
charge. Other studies deal with nuclear underground 
explosions and analyse the records in terms of their 
power spectral density [55]. 

The air blast effects are mainly related to the ex- 
plosion of nuclear weapons in the atmosphere. The pres- 
sure wave resulting from an explosion near ground con- 
sists of an abrupt rise in pressure followed by a decay 
from which a negative pressure half-wave results [56]. 
The shape of the pressure wave is well defined and so 
its effects on structures can be studied by using the 
deterministic theory of vibrations. 

The value of the load to be adopted in design has 
to be established on a strategic basis taking into account 
the degree of protection that is desired. 

Sonic boom can also be considered as a type of 
blast load. 


5 — CONCLUSIONS 


The main purpose of this introductory report is to 
serve as a basis for the discussion to be held during 
the Congress. So the present conclusions contain propo- 
sals on research subjects about which discussion is 
desirable. 


5.1 — The need of a correct definition of the loads 
acting on the structures was emphasized. It is important 
to state this definition in a basic way according to well- 
established general theories. Discussion on the most 
convenient methodology to attain this scope is welcome. 


5.2 — Recent progress concerning the knowledge 
of the dynamic action of wind was described. 

To increase the available information it seems 
desirable: 


5.2.1 — To record systematically, by convenient 
anemometers, the turbulence of wind, in supplement to 
the actual recording of mean velocities. 


5.2.2 — To improve the actual representation of 
wind by using velocity spectra that duly vary in function 
of geographical conditions and other pertinent variables. 

5.2.3 — To establish on sound statistical basis the 
velocities to be used in design. 

5.2.4 — To study the dynamic behaviour of structu- 
res taking simultaneously in consideration the turbulence 
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inherent to the natural wind and the turbulence created 
by the structure itself. For this purpose model studies 
considering the wind velocity profile and turbulence 
effects seem promising. 


5.2.5 — To pursue the observation of important 
structures in order to get more information concerning 
their behaviour under wind actions. 


5.2.6 — To include in the building codes not only 
simplified data for the design of usual structures, but 
also basic data that can serve for desing of important 
structures. 


5.3 — For further progress in the definition of seis- 
mic loads the following lines seem promising: 


5.3.1 — To continue installing strong-motion accel- 
erographs and duly interpret the obtained records. 


5.3.2 — To define the seismic loads by the spectral 
density of acceleration and consider the variations of the 
spectra due to local conditions. 


5.3.3 — To define seismicity of a region by the 
probability of a given level of spectral density of accel- 
eration being reached in that region. The seismicity of 
the different regions must be established by combining 
geophysical, geological and seismic information and 
using, as far as possible, quantitative statistical criteria. 


5.3.4 — To pursue analytical and experimental stu- 
dies, based on the theory of random vibrations, mainly 
for obtaining further information concerning the beha- 
viour of linear and non-linear systems with several de- 
grees of freedom. 


5.3.5 — To complement the installation of strong- 
motion accelerographs for recording soil accelerations, 
by installing equipment also allowing to observe the 
behaviour of the structures themselves. Useful informa- 
tion may also be obtained by dynamic tests of real 
structures, even under vibrations of small amplitude. 


5.3.6 — The improvement of building codes is most 
desirable. As for wind, it would be convenient that codes 
contain simplified rules for the design of ordinary struc- 
tures and basic data to be used in special studies of 
important structures. 


5.4 — Discussion is also open on dynamic loads 
other than wind and earthquakes. Among these, loads 
due to traffic, machinery and blast were mentio- 
ned. New data established on modern scientific bases 
shall largely contribute to the design of more economical 
and safer structures. 


5.5 — Finally, the importance of international col- 
laboration as a powerful means of accelerating progress 
must be emphasized. This collaboration may be partic- 
ularly fruitful for the establishement of recommenda- 
tions of general character on which regional codes may 
be based. 
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SIGNIFICADO DO CAMPO ELÉCTRICO NO ARVOREJAMENTO 
DE POLÍMEROS DEVIDO A UMA HETEROGENEIDADE 
ELÍPTICA E HIPERBÓLICA 


RESUMO 


Calcula-se a distribuição espacial do campo eléc- 
trico em torno de um elipsóide e de um hiperbolóide de 
revolução artravés das linhas (superfícies) equipoten- 
ciais, que definem a direcção do campo eléctrico, e in- 
troduz-se a noção de linhas (superfícies) isocampo para 
caracterizar a intensidade do campo eléctrico. A compa- 
ração da teoria com a experimentação do arvorejamento 
de polietileno permite interpretar o fenómeno pela acção 
do campo eléctrico. 


LISTA DE SIMBOLOS 


A — área de uma secção; 

a — semi-eixo maior de um elipsóide (ou elipse) 
e parâmetro de um hiperbolóide (ou hipér- 
bole); 

a; — semi-eixo maior do eléctrodo elíptico; 

b — semi-eixo menor de um elipsóide (ou elipse) 
e parâmetro de um hiperbolóide (ou hipér- 
bole); 

b E — semi-eixo meno: do eléctrodo elíptico; 


by — parâmetro do eléctrodo hiperbólico; 


c — distância focal comum aos elipsóides (elip- 
ses) e hiperbolóides (hipérboles); 


d — espessura do dieléctrico; 
E — intensidade do campo eléctrico; 
E — campo eléctrico; 
Ec max — intensidade máxima do campo eléctrico no 
caso de um eléctrodo elíptico; 
E 4 mas — intensidade máxima do campo eléctrico no 
caso de um eléctrodo hiperbólico; 
Eo — intensidade do campo eléctrico uniforme; 
Entao intensidade máxima do campo eléctrico; 


(*) Recebido na redacção em 14-9-971 
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SYNOPSIS 


We calculate the spatial distribution of electrical 
field with an ellipsoidal and hyperbolical electrode by 
means of equipotencial lines (or surfaces), which define 
tha field direction, and the new notion of isofield lines 
(or surfaces) is introduced to describe the field intensity. 
A comparison of this lines with the treeing of polyethy- 
lene conduct us to an interpretation of the phenomenon 
by the electrical field action. 


E (z) — intensidade do campo eléctrico no eixo do 
eléctrodo heterogéneo; 


=p —p 


E,EG Es — componentes do campo eléctrico nas direc- 
dos eixos coordenados nr, p ev; 
9, 9, 9, — coeficientes métricos; 
Kokim kom k pkq — constantes; 
| — comprimento; 


Pin), Q,(7) — soluções da equação de Legendre de 1. 
ordem; 
r — raio de curvatura da ponta do eléctrodo he- 
terogéneo; 
re — raio de curvatura da ponta elíptica; 
rw — raio de curvatura da ponta hiperbólica; 
s — distância disruptiva; 
V — tensão (valor eficaz); 
x,y, Z — coordenadas carteseanas; 
x,y, z — vectores unitários do sistema carteseano; 
x — parâmetros das rectas da zona de transição 


da função E ., (s); 
v — função potencial; 
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Yp (4) — função potencial particular; 


7%, 0, “ — coordenadas do sistema elíptico; 
ego cmi cmi 
7%, 0, > — vectores unitários do sistema elíptico de 
coordenadas; 
re — parâmetros do eléctrodo elíptico; 
0 | — parâmetro de um hiperbolóide (hipérbole); 
1y — parâmetro do eléctrodo hiperbólico; 


: — distância de um ponto do espaço ao eixo z; 


+ — volume. 


1. INTRODUÇÃO 


A crescente aplicação de materiais poliméricos na 
electrotecnia e em especial de polietileno no isolamento 
de cabos eléctricos para altas frequências e de altas 
tensões mostraram significativa importância no chamado 
teste de agulhas [1], dada a grande sensibilidade do 
polietileno a fenómenos de ionização e porque nestes 
ensaios laboratoriais a energia se pode controlar cômo- 
damente e assim evidenciar a fase anterior à disrupção 
completa do material, o que não é fácil sob condições 
normais de serviço. 

De facto o fenômeno de arvorejamento (do inglês 
«treeing») tem sido sempre analisado, teórica e experi- 
mentalmente, através de um sistema hiperbolóide-plano 
devido à realização dos ensaios de investigação 
se simplificar pela introdução de agulhas de coser em 
aço niquelado no dieléctrico. Um sistema elipsóide-plano 
representa no entanto em frequentes circunstâncias um 
modelo de análise mais rigoroso, como sucede quando 
se aplica uma tensão entre dois eléctrodos paralelos 
com uma heterogeneidade adjacente a um deles e de 
comprimento muito menor que a espessura do dieléc- 
trico. É exemplo evidente um cabo eléctrico em que, por 
defeito de fabrico, apresenta uma incrustação condutora 
no material isolante, quer do lado do condutor quer do 
lado da armadura. 

O presente trabalho ocupa-se esencialmente do cál- 
culo das distribuições do campo eléctrico em eléctrodos 
com perfil elíptico e hiperbólico e da comparação com 
o desenvolvimento do processo de arvorejamento. Além 
de se efectuar o cálculo em torno de um elipsóide metá- 
lico analisa-se também a distribuição do campo eléctrico 
à volta de um hiperbolóide (na literatura encontra-se o 
cálculo apenas no eixo da respectiva hipérbole, como 
formulou Mason em [2]). 


2. HETEROGENEIDADE ELIPTICA 


Considera-se fundamentalmental o estudo da distri- 
buição do campo eléctrico no espaço que rodeia uma 
configuração elíptica. Para se alcançar este objectivo 
admite-se isotropia do dieléctico e uma configuração 
geométrica ideal dos eléctrodos. 


2.1. Distribuição do potencial 


De um modo genérico a função « (x, y, Z) represen- 
tativa da distribuição do potencial num espaço caracte- 
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rizado pelo sistema de eixos carteseanos (x, y, Z,), satis- 
faz à equação de Laplace 


+ “ s 
dºw d*w dy 


á Nom N 5 (1) 
dx” dY- dz” 
Se nesse espaço reinar um campo elétrico uniforme 
Eo. com uma dada direcção, as superfícies equipoten- 
ciais são planos normais às linhas de campo (fig. 1 a). 
Se se introduzir porém um volume metálico no espaço 
considerado a distribuição do potencial é alterada em 
iunção da forma dessa heterogeneidade, de modo que a 
superfície do volume metálico seja uma superfície equi- 
potencial e a uma distância infinita a perturbação seja 
nula. Assumindo que o volume metálico tem a forma 
geométrica de um elipsóide de revolução, com o eixo 
maior segundo o eixo de coordenadas z e o centro na 
origem de coordenadas, o plano x, y, normal ao elipsóide 
e passando exactamente pelo centro, é também uma su- 
perfície equipotencial (fig.1 b). 


q) 


c) d) 


Fig. 1 — Esquematização do princípio de cálculo de 
uma configuração elíptica (linhas equipo- 
tenciais a traço-traço e linhas de campo 
eléctrico a traço contínuo). 

a) Campo uniforme. 

b) Perturbação do campo devida a um 
elipsóide metálico. 

c) Formação do eléctrodo elíptico por um 
plano e um semi-elipsóide. 

d) Configuração elipsóide-plano. 
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GESTRA e GESTRA-KSB 


Purgadores de funcionamento termostático 
e termodinâmico, tipos BK, HK e UNA- 
MAT PN 6/10/16/25. 

Construção normal: corpo de aço forjado 
C 22, órgãos interiores de aço inoxidável. 
Trabalham automáticamente até 8/17/22 
atmosferas de pressão diferencial. Liga- 
ções soldadas, roscadas segundo BSP ou 
NPT e com flanges segundo DIN ou ASA. 


BK 15 HK 
UNAMAT 


Execuções para altas 
mesmo principio de 


pressões 
funcionamento até 
PN 500 (existentes em PN 40/100/160 


com o 


/320/500). 


BK 17/18/110/111/112 


Purgadores termodinâmicos para grandes 
caudais até 200m'/h e para pressões 
máximas, com órgãos de obturação esca- 
lonadas GESTRA, até PN 640. 


GK/NK 


a ans 


L] 


Válvula de descarga GESTRA-KSB tipo AK. 
A válvula de descarga tipo AK serve 
para a purga no arranque de grandes ins- 
talações e turbinas. Fecha-se ao atingir 
certa temperatura sem energia auxiliar. 


Construções até PN 500 (existentes em 
PN 40/160/500). 


AK 


[) 


Purgadores de flutuador GESTRA-KSB 
com dispositivo para purga de ar auto- 
mático ou manual, para qualquer apli- 
cação até PN 160 (existentes em PN 
16/25/40/100/160). Construções espe- 
ciais para a purga de redes de ar compri- 
mido. 


UNA-Spezial 


Válvulas de retenção GESTRA, de Disco, 
tipo RK, nos diâmetros nominais de 
15-400 mm até PN 40. Corpos de 
bronze, aço vazado, aço inoxidável, etc. 
Órgãos interiores sempre de aço inoxi- 
dável. Para fluidos voláteis, com disco 
revestido a Viton ou Perbunan. Às vál- 
vulas CESTRA de disco tipo RK impõem- 
-se sobretudo pelas suas reduzidas dimen- 
sões, pela sua vedação perfeita, peso 
minimo e fácil montagem. 

RK 46 


AK 41/44/46 


[| 


Visores VKN e VK 35 para o controle vi- 
sual da passagem de todos os tipos de flui- 
dos em tubagens para o controle de per- 
mutadores de calor e do funcionamento de 
purgadores (existentes em PN 16/25/40). 


VKN 


[) 


Válvulas misturadoras GESTRA, Calormix 
MT 22, para a produção de água quente 
por mistura de água fria com vapor. Regu- 
lação manual sem graduação. Ruídos mi- 
nimos. Ausência de resíduos de calcário 
mesmo a temperaturas críticas, por distri- 
buição adequada da câmara de mistura. 


MT 


L) 


Reguladores automáticos de temperatura 
de saída, tipo BW, com comando termos- 
tático para adaptação a instalações de 
aquecimento, com água quente e sobrea- 
quecida e com termofluidos. 

Válvulas limitadoras de temperatura de 
saída de água de refrigeração, tipo CW, 
com comando termostático para a redução 
do consumo de água de refrigeração. 


Válvulas de fecho rápido para a purga 
periódica de caldeiras, tipos ZA, SPA, PA, 
com accionamento manual ou por pedal. 
Aptas para todas as condições de serviço 
e para todos os tipos de caldeiras. Estas 
válvulas existem também com comando 
a distância e programação eléctrica (exis- 
tentes em PN 25/40/64/160/250). 


Válvulas reguladoras GESTRA Reactomat 
BA, com órgãos de obturação escalonados 
para a purga continua de lamas em cal- 
deiras (existentes em PN 40/64/100/ 
/160/250/320). 


BA 


[) 


Filtros SZ até DN 50 (27) e PN 500 
(existentes em PN 40/160/320/500). 


LU) 
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Sob o ponto de vista electrostático a distribuição 
de potencial não se altera se se substituir uma superfície 
equipotencial por uma superfície metálica. Deste modo 
a distribuição de potencial considerada mantém-se se se 
considerarem dois eléctrodos constituídos por planos ili- 
mitados paralelos e infinitamente afastados um dos quais 
contém um semi-elipsóide com o eixo maior normal ao 
respectivo plano exactamente na origem de coordenadas 
admitidas (fig.1 c). 

Assim uma tensão V aplicada entre os eléctrodos 
estabelece um campo eléctrico no espaço compreendido 
entre eles constituído pela sobreposição de um campo 
uniforme E, e um campo de perturbação introduzido pelo 
elipsóide (base do pensamento da solução (20 a) adian- 
te proposta). 


O estudo da configuração geométrica estabelecida 
é simplificado num sistema de eixos coordenados trior- 
toganais (7, 5, vY ) caracterizados pelos elipsóides de 
revolução de parâmetro 7 e confocais com o elipsóide 
ng que forma a heterogeneidade condutora ( np <n<. 00) 


z” 2 


SE a (2) 


LR n— 1 


onde : representa a distância de um ponto do elipsóide 
ao eixo z e c a distância focal comum (fig. 2 a); pela 
família de hiperbolóides de parâmetro à normais aos eli- 


psóides anteriores (0< < 7) 
zº p? E 
E é PR — €? 3 
6? 1—6 


e pelos planos contendo o eixo z e fazendo o ângulo “ 
(em que O <v<'7) com uma dada referência, por exem- 
plo o plano x, y. 

Um ponto definido por este sistema de coordenadas 
(n, 1, Y) tem uma distância ao eixo z dada por 


e-eV (RN) (16) (4) 


como se deduz de (2) e (3). 
Dada a simetria das figuras de revolução em torno 
do eixo z basta estudar as elipses no plano x, z 


1 (5) 


com os semi-eixos a e b tais que 
c'—- a — b? (6) 


caracterizando-se cada eclipse pela excentricidade 


pin (7) 
a 


Nestas condições as hipérboles normais respectivas 
são caracterizadas por q satisfazendo a importante pro- 
priedade 


Zz = cr6 (8) 
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derivada de (2) e (3), onde a distância focal comum 
c é estabelecida pelo eléctrodo elíptico através dos seus 
semi-eixos maior a ; e menor b ; 


e-Vaç—b, (9) 


conforme com (6). 

Para se exprimir a equação de Laplace (1) no novo 
sistema de eixos coordenados considere-se um compri- 
mento elementar do espaço carteseano 

(dlZ= (dx + (dy) +-(dz) (10 a) 


que se transforma para o novo sistema de coordenadas, 
com a condição necessária e suficiente de ortogonalidade 


(ai) = oi (60) + as(a0) + ai (02) 
onde os coeficientes métricos são 
pi | 3 ES . = k 


e = ) + + 


(10 b) 


ev tda) E ] 


A partir das relações entre as coordenadas dos dois 
sistemas em estudo 


x = c shyr sinf cosy 
y = eshr sinô siny (12) 
Zz = c chm cosº 
obtém-se então 
Im ini Cc / SÃ E d 
V a 


(13) 


— p a 


Com os eixos unitários carteseanos x, Y, Z pode 

escrever-se a definição 
do — 

grad q = ds E qu = 

dx dy dz 


z (14 a) 


- jp ep 
e com os vectores unitários q 0 “Y do sistema de 


coordenadas elíptico 


dia ps! RA A. + ' BE dim 
Gm dn 9 db 04 dy 
(14 b) 
atendendo a (11). 
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Pela definição de divergência de um vector 


div E — lim Ee (15 a) 


sendo o volume 47 determinado pelos comprimentos 
0; 47, 90 4h e gy; ay (íig. 2 b) 


à7 — 9, 96 9 dn Al AS, (16) 


Nestas condições a divergência bruta do campo 
eléctrico em relação à superfície A, -— 9541 9, A, e 


—+ 
segundo a direcção xr tem a expressão 


RU 


9, 9, E, 


E ei , 
| 9, 8, E, d dt (9, 9, E, ) ax 


Ar, 
Ir, 


E) + À 
| (9, 9, E. ) MO A 
f | o) 


e análogamente para as restantes superfícies Ay e A; do 
volume Ar, pelo que 


— —» d | » N 7 e À 
PERA | dn (9,9,E,) df (9,9, E, 


de SO A 


Ea E 
dy (9, 9 ? 


e, através de (16), 


(15 b) 


Deste modo o laplaceano da função potencial, se- 
gundo a definição 


Ay div grad « (17 a) 


resulta 


, 1 d [90 94 E Ô 
Ag .— o 
8,94 8, dr 9, dr 


o (+ 4 d [9a 9 e) 
db In 09, dy 9, d + 


(17 b) 


Estas expressões são gerais, válidas para qualquer 
sistema de coordenadas triortogonais. 


No sitema elíptico definido, por razões de simetria, 
as grandezas eléctricas são independentes de >, isto é, 
de | dy O |, resultando portanto a equação de Laplace 


» [9994 dy ) [n9, 
à (0 dm), à (mo deb q 
dm 9, dr d9 9, dh 
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Fig. 2 — Sistemas de coordenadas das configurações. 
a) Sistema elíptico de coordenadas. 


b) Volume para o cálculo de div E. 

c) Sistema de estudo da configuração 
elíptica. 

d) Sistema de estudo da configuração 
hiperbólica. 

e) Configuração elíptica experimental. 

f) Configuração hiperbólica experimental. 


dando as expressões (13) finalmente 


oe LN Fon. do (A 


ALR 
equação diferencial que se integra por reflexão física 
introduzindo uma solução co tipo 


" Eo 2z | kúcp (x) (20 a) 


once a primeira parcela atende ao potencial correspo-- 
cente ao campo uniforme E; e a segunda parcela consi- 
dera a não uniformidade consequência da heterogenrei- 
dade elíptica, sendo k uma constante arbitrária e «<p (x) 
uma função particular apenas de 7 a determinar. 


Por substituição de (20 a) em (19) obtém-se 


d* vo (7) dp (rf) 
dá p 
O — 3) ———— = 2if= 


dn” 


+ 2vp (7) Õ 
(21) 


dt 


que é uma equação de Legendre de 1.º ordem. Como 
7) > 1 a solução desta equação é [3] 


vp (1) kp Pilr) + ko O (m) 
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onde kp e ka são constantes arbitrárias e 


P(n) = 2 
% 7 + 1 
Q, (7) do, E Ee 
dis 2 q — 1 
ou seja 
ken k AR a CL 
Cplr) CA nesta Mas - k 
ni F a ( 2 e À ) 


A distâncias infinitas o potencial deve ser constante 
e independente da perturbação introduzida pelo elipsóide, 
pelo que 


lim vplr) — o 
%—+ OC 


condição que permite afirmar kp = O, e portanto com 
(20 a) 


Es e dt tef lã is 1) (20 b) 
= Z — e n e 
' ” O cç 2 n— 1 
e supondo nulo o potencial do eléctrodo elíptico, o que 
em nada restringe o raciocínio, determina-se 


Eq € 
k ka é DE: CR 
1 te +1 1 
In — - 
2 te — 1 rE 


ea EA ads 
é ig] Hom eim (22) 
1 TE +41 1 
, 2 ve “Ve 


2.2. — Distribuição do campo eléctrico 


O campo electrostático calcula-se pela lei funda- 
mental 


—+, 


E grad « (23) 
e através de (22), (14 b) e (13) 


1 % r 
; Eme ) 
fe — 4 : 2 n— 1 n—1 
EDS fi À eme ue r Y 
Eo y KR — (1º kr *E 41 x 1 
N 2 E —1 TE 
Att (24) 
Ly poa 1 E - A > 7] 
| 2 — à? 1 ZE +-1 1 


ms fi 
2 te —1 Lá = 

Sobre o eixo z (quando 6 = 1 o hiperbolóide dege- 

nera na recta correspondente ao eixo Z) a componente do 


campo eléctrico segundo 1 anula-se, donde 


/ 1 r 1 7 
sia in — 
E(z) 2 q = 1 7º 1 —», 
1 z 
Eo 1 HE+-1 1 
| Treat ——— 
2 te — 1 TE 


(25) 
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tendo o campo eléctrico a máxima intensidade na ponta 
do elipsóide relativo ao eléctrodo (para n = me ) 


1 “E 4-1 "E 
SE a E 
Emax 2 “E —1 ng *—1 
EEE + (26) 
Eo 1 te +1 1 
esmo * Rim es 
2 TE -— nE 


Dividindo (25) por (26) obtém-se a distribuição 
do campo eléctrico no eixo z normada em relação 
à Eax* OU seja 


da dará fat ei 
Caso PR -E | 2 €-n+T 
1 n 
o ne há $ E À . (27) 


Para grandes distâncias da ponta do eléctrodo 
elíptico 


lim E(z) = Es (28) 


Z —> 00 


como faácilmente se vê de (25), tendendo-se rápidamente 
para a uniformidade (fig.3 a). 

Observa-se que os cálculos efectuados supõem a 
permitividade do meio ambiente constante, desprezan- 
do-se portanto a influência de efeitos de polarização 
dieléctica, e além disso são válidos para dois eléctrodos 
infinitamente afastados. Relativamente à primeira simpli- 
ficação acrescenta-se que no caso de dielécticos sólidos 
é em geral justificada, sobretudo em materiais não pola- 
res, como é o caso do polietileno. Em relação à segunda 
hipótese mostra a expressão (22) deduzida que as equi- 
potenciais são práticamente paralelas para distâncias da 
ponta do semi-elipsóide da ordem de 3a, (fig. 3 a), ou 
seja, considerando um eléctrodo plano afastado de 
d=4a- do outro que contém o semi-elipsóide (fig. 
1 d) a distribuição do campo eléctrico não é prática- 
mente perturbada pois a intensidade do campo difere 
do valor Ey a menos de 1% (fig. 3). 

Nos casos em que a distância disruptiva entre a 
ponta do eléctrodo heterogéneo e o eléctrodo plano 
oposto 


s=d-— ar (29) 


é muito pequena esta segunda aproximação deixa de ter 
significado, porquanto as superfícies equipotenciais se 
afastam grandemente da planaridade junto do eixo do 
eléctrodo elíptico, devendo considerar-se outro critério 
ce cálculo. Quando a ponta do eléctrodo de revolução 
está muito próxima do eléctrodo oposto (condição 
s << ap) o comportamento da configuração afigura-se 
porém como o de um hiperbolóide. 


3. HETEROGENEIDADE HIPERBÓLICA 


Abstraindo do estudo particular de um eléctrodo 
hiperbólico somos conduzidos à analise do modelo hiper- 
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Fig. 3 — Configuração elíptica: 
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r;=0,1 mm,aç = 0,9 mm. 


a) Linhas equipotenciais, segundo a eq. (22); parâmetro normado «q / Es. 
b) Linhas isocampo, segundo a eq. (24); parâmetro normado E/E,. 


bolóide - plano para interpretar os casos em que a dis- 
tância disruptiva de uma heterogeneidade elíptica é mui- 
to reduzida. 


3.1. — Distribuição do potencial 


Considere-se um sistema constituído por um eléc- 
trodo com a forma de um hiperbolóide de revolução com 
o eixo normal a um outro eléctrodo plano (fig. 2 d). 

Seja c a distância focal do hiperbolóide definido 
num sistema casteseano de coordenadas pelos valores 


seb, tais que 
c = V sº + by (30) 


onde s é o afastamento da ponta do eléctrodo até à 
origem, podendo caracterizar-se o eléctrodo pelo parã- 
metro 


b 
'yu= arc tg a (31) 


Aplicando uma tensão V ao hiperbolóide define-se 
uma distribuição do campo eléctrico no dieléctrico de tal 
maneira que as linhas de campo sejam normais às su- 
perfícies de ambos os eléctrodos. 
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Dada a configuração geométrica do sistema pre- 
fere-se também efectuar o seu estudo num sistema de 
coordenadas ( 7, 1, « ) definido pelos hiperbolóides 
9(fy< 0<. =) 2) confocais com o eléctrodo hiperbó- 
lico, os elipsóides 7 (0<.r<.oo ) normais aos hiperbo- 
lóides anteriores e os planos v(0<“U< =) contendo 
o eixo z e fazendo o ângulo “ por exemplo com o plano 
x, z. Observada a simetria em relação ao eixo z basta 
analisar o modelo hiperbolóide - plano no plano x, z, isto 
é, independentemente de “, reduzindo-se as figuras de 
revolução a hipérboles e elipses. 


Assim um ponto de coordenadas carteseanas (x, Z), 


existente sobreuma hipérbole de parâmetro 9 definido 
poraeb 


0 = arc tg , (32) 
satisfaz a equação da hipébole 
= E 1 (33) 
a? b? , 


Fisicamente as linhas equipotenciais da configura- 
ção hiperbólica são hipérboles, donde se conclui a in- 
dependência da função potencial em relação a 7 (além 
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de não depender “, como se referiu). Nestas condi- 
ções a equação de Laplace (17) simplifica-se 


od [a 9 d 
(E. l-o (34) 
dó 96 do 
e por intermédio de (11) e (12) obtém-se 
ds d: 
E + te = 6. (35) 
de” 
cuja solução [4] é de forma 
Y = k,H + KH In cot (36) 


onde as constantes k, e k | se calculam pelas condi- 
ções limites 
1) 8 = 


0H , 2 Wes 


2) 0 = - 14 = 0; eléctrodo plano à terra; 


N2z 


AVI 
d 


Ny 


| adido, 


| 

má ANG sauna ne “UM. au 

| SA CARL RENM 
NGN AEE a 
+ Rm 

E 

MIS 

K mm Enma nes lista) 
EE 


Fita fondo mto 


a 


ALL LLA 


resultando finalmente 


(37) 


As hipérboles relativas às equipotenciais serão ca- 
racterizadas por conseguinte pelo parâmetro 


a(s" 


uma vez normado o potencial com a tensão V apli- 
cada entre os eléctrodos (fig. 4 a, à esquerda). 


*w (38) 


2 


1 


ae 


x [mm] — o 


1 .. 
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Fig. 4 — Comparação das distribuições do campo eléctrico nas configurações elíptica (em cada semi- 


-gráfico direito) e hiperbólica (à esquerda): d = 


4,5 mm, s=3 mm, r= 0,01 mm. 


a) Linhas equipotenciais (traço contínuo) nmormadas com a tensão V e linhas de campo estimadas 


(traço-traço). 


b) Linhas isocampo normadas com E ax 


c) Intensidade do campo eléctrico no eixo do eléctrodo heterogéneo. 
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3.2. — Distribuição do campo eléctrico 


Seguindo igual raciocínio a equação (22) permite 
calcular a intensidade do campo eléctrico 


V 
E Ee 
Cy (39) 
c sin 6 | sh? r + sin'6 In cot 
onde 
: b 
sin 1 nec (40) 
c 
como se vê facilmente (fig. 2 d), e 
h (41) 
ug 
b 


deduzida da primeira equação (12) e de (40). 
O valor máximo do campo eléctrico verifica-se na 
ponta do eléctrodo hiperbólico (para 1 ty er =0) 


V 


max ns 'y (42) 
c sin* 6, - In cot 


pelo que a expressão normada da intensidade do cam- 
po resulta 


E sin” ty 
a (43 a) 
E max sin 0 ysh"n + sin” 9 
ou, com (40) e (41), 
E b'W 
RT Sá (43 b) 
E max Vb' + c* x? 


Esta expressão permite calcular o valor da abcissa x 
de um ponto sobre a hipérbole caracterizada por b (ou 
seja 1 ) onde reina um dado valor E/E ,,.,: com a equa- 
ção (33) calcula-se a correspondente ordenada z (fig. 
4 b, à esquerda). 

A distribuição da intensidade do campo eléctrico 
go longo do eixo do eléctrodo hiperbólico obtém-se de 


(39) fazendo «n=o 
V 
E (z) 5 Zá 
c sin" 6. In cot dé.) a 
2 
donde 
E(z) sin” 6y 
= — 45 
E max sin” 6 a) 
ou 
E(z) b'y 
= = . (45 b) 
Emas b* 


É importante notar que estas expressões deduzidas 
para o potencial e para o campo eléctrico são indepen- 
dentes da espesura do dieléctrico, isto é, admitiu-se que 
todo o espaço é ocupado pelo mesmo meio ambiente. 
No caso real em que o dieléctrico onde se encontra a 
ponta do hiperbolóide é limitado haverá uma diferença 
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por efeito da superfície da transição com o ar, que será 
tanto menor quanto mais reduzida for a permitividade 
do material isolante. 

As grandezas definidas para o campo eléctrico po- 
cem normar-se facilmente por intermédio da intensi- 
dade do campo eléctrico uniforme E, que existiria se 
embos os eléctrodos fossem planos e tendo em atenção 
a espessura d do dieléctrico 


E ki 46) 
o 3 ( 
mas preferiu-se a normação com E... por Ego não ter 
significado real num sistema hiperbolóide - plano. Este 
pensamento justifica a mesma normação no caso do 
modelo elipsóide - plano a fim de se efectuar uma com- 
paração entre as duas configurações. 


4. COMPARAÇÃO DAS HETEROGENEIDADES ELÍPTI- 
CA E HIPERBÓLICA 


O cálculo numérico das expressões deduzidas, rea- 
lizado num computador digital, permitiu uma represen- 
tação gráfica das grandezas eléctricas com a vantagem 
de uma visualização clara das diversas distribuições. 

O tratamento analítico da configuração elíptica foi 
realizado com a origem de coordenadas no centro do 
correspondente elipsóide de revolução (fig. 2 c), de mo- 
do que o campo eléctrico decresce a partir da ponta do 
eléctrodo à medida que os valores da ordenada z aumen- 
tam; na configuração hiperbólica porém o eléctrodo en- 
contra-se afastado da origem (fig.2 d), decrescendo a 
intensidade do campo desde a sua ponta e à medida que 
os valores de z diminuem. Para comparar as duas confi- 
gurações no mesmo sistema de coordenadas carteseanas 
e para iguais distâncias s do eléctrodo plano até à pon- 
ta do eléctrodo heterogéneo exige-se uma mudança de 
coordenadas. Tomando os eléctrodos heterogéneos na 
posição de ánodo, isto é, com os correspondentes eléc- 
trodos planos ao potencial nulo (fig. 2 e, f), os cálculos 
da configuração elíptica foram convertidos no sistema 
de coordenadas do modelo hiperbólico tendo em consi- 
deração a espessura d do dieléctrico. 

Quando os eléctrodos elíptico e hiperbólico pos- 
suem uma ponta com igual raio de curvatura, respecti- 
vamente 


r nd 
p= - 47) 
dg ( 
LS) 
bº, 
FH = (48) 
s 


o último contorno define um eléctrodo ligeiramente mais 
grosso (comparem-se os respectivos perfis na fig. 4). 
Mas pode assumir-se que o primeiro eléctrodo corres- 
ponde concepcionalmente ao segundo no qual se rodou a 
superfície do hiperbolóide de maneira a formar o plano 
onde assenta o semi-elipsóide da configuração elíptica. 
As diferenças introduzidas na distribuição do campo 
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eléctrico são devidas essencialmente à existência deste 
plano (que não existe na configuração hiperbólica). 


1. A direcção do campo eléctrico traduz-se pelas 
equipotenciais (fig. 4 a), que lhe são constantemente 
normais. Nos dois casos estudados constata-se uma dife- 
rença na direcção do campo eléctrico nos pontos de igual 
abscissa sobre os planos transversais ao eixo do eléctro- 
do heterogéneo, sobretudo junto da sua ponta: a compo- 
nente longitudinal é maior que a transversal em ambas 
as configurações mas a relação entre estas componentes 
é menor no caso da ponta hiperbólica, quer dizer, a aber- 
tura das linhas de campo é maior se se trata de um 
eléctrodo hiperbólico. 


2. As linhas isocampo, definidas como sendo uv 
lugar geométrico dos pontos com igual intensidade de 
campo eléctrico (correspondendo no espaço a superfícies 
isocampo), caracterizam a intensidade do campo eléc- 
trico (fig. 4 b). Em primeiro lugar o seu desenvolvimento 
ao longo do eixo do eléctrodo heterogéneo é muito mais 
comprido para a configuração elíptica; em segundo lugar 
o desenvolvimento em planos transversais a esse eixo na 
zona de disrupção é mais largo para o eléctrodo elíptico 
(embora em planos que atravessem o eléctrodo possa 
ser menor). 


3. Para uma dada tensão a intensidade máxima do 
campo eléctrico é maior na ponta hiperbólica do que na 
elíptica: se r= 0,01 mm, s=3 mm e d=4,5 mm 
obtém-se os valores E | max/Eo = 67 a partir de (26) 
e Ey max/Fo = 127 através de (42), portanto numa re- 
lação de 1:2 (fig. 4 c, para z = 3 mm). Outro modo de 
formular este facto: para iguais valores máximos da in- 
tensidade do campo eléctrico a tensão aplicada ao eléc- 
trodo elíptico deverá ser inferior ao correspondente do 
eléctrodo hiperbólico (cerca de metade no exemplo cal- 
culado). 

4. A intensidade do campo eléctrico sobre o eixo 
do electrodo (fig. 4 c) reduz-se muito rãpidamente à 
medida que se afasta da ponta em ambos os casos, mas 
enquanto na configuração hiperbólica tende para um va- 
lor muito reduzido no eléctrodo elíptico tende para o 
valor Eo do campo uniforme. 

5. A distância disruptiva, ou seja, o afastamento 
entre a ponta do eléctrodo heterogéneo e o eléctrodo pla- 
no oposto, é decisiva na determinação do valor máximo 
da intensidade do campo eléctrico (fig. 5). Numa confi- 
guração hiperbólica se s é muito pequena o campo na 
ponta tem um valor muito elevado; para maiores distân- 
cias s o campo máximo reduz-se tendendo para um an- 
damento pouco decrescente quando a distância disrupti- 
va aumenta. Por seu lado numa configuração elíptica 
admitiu-se que o eléctrodo oposto ao eléctrodo elíptico 
se encontrava infinitamente afastado pelo que para pe- 
quenas distâncias s os cálculos deixam de ter signifi- 
cado; substitui-se então o raciocínio de cálculo pela con- 
figuração hiperbólica equivalente (mesma distância dis- 
ruptiva e ponta com igual raio de curvatura), Segundo 
este critério há uma zona de transição que se define 
numa representação gráfica semi-logarítmica pelas rectas 
tangentes às curvas relativas á função E ax (s) (ou à 
função normada E, .., (s)/V, pois se facilita o cálculo 
para diferentes tensões aplicadas) com igual parâmetro 


TÉCNICA 413 


r. No caso em que d = 4,5 mm (fig. 5) é válido o cál- 
culo da configuração elíptica até s = 2 mm; abaixo de 
s = 0,75 mm reinam os cálculos da configuração hiper- 
bólica equivalente; entre s=2 mm e s= 0,75 mm a 
zona de transição é caracterizada pela recta 


E 
log mondo a — 0,435 (49) 
V 
ou 
Eai V; 10 0988 (50) 
com E ax expresso em MV cm-le V em KV, sendo « 


uma constante dependente do raio de curvatura, como 
a seguir se indica para o exemplo calculado: 


r x 
(mm) (MV em! /KV) 
0,01 0,45 
0,02 0,26 
0,10 0,08 


6. O raio de curvatura da ponta do eléctrodo hete- 
rogéneo é igualmente de grande importância na definição 
da intensidade do campo eléctrico e em especial do valor 
máximo (fig. 5). Quanto menor for r maior é a inten- 


Fig. 5 — Intensidade máxima do campo eléctrico 
na ponta do eléctrodo heterogéneo em fu- 
ção da distância disruptiva s (coordena- 
da da ponta); d = 4,5 mm. 
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sidade para igual distância disruptiva, sendo o corres- 
pondente valor da configuração hiperbólica sempre su- 
perior ao da elíptica. 


5. INTERPRETAÇÃO DO ARVOREJAMENTO PELO 
CAMPO ELÉCTRICO 


O estudo dos modelos definidos tem importância 
por possibilitar um esclarecimento físico do arvorejamen- 
to de dieléctricos sólidos, conforme a experimentação 
revela. 

O fenómeno de arvorejamento é observável ôptica- 
mente nos materiais isolantes poliméricos, particular- 
mente nos que possuem elevada transparência, como por 
exemplo o polietileno. Este material foi utilizado na in- 
vestigação experimental, para o que se embeberam eléc- 
trodos metálicos com os perfis examinados, imergin- 
do-se o conjunto num banho de óleo de transformador 
para evitar descargas superficiais à tensão do ensaio. 

Quando a tensão aplicada ao dieléctrico for sufi- 
ciente para provocar localmente, em virtude da presença 
da heterogeneidade metálica, uma intensidade de campo 
eléctrico que origine a disrupção intrínseca incompleta 
do material forma-se uma cavidade onde surgem des- 
cargas parciais e daí o desenvolvimento de canais que 


a) b) 


t aé 06 80 Mt É 


mm e jt | 


correspondem ao arvorejamento. Estes canais, sendo 
devidos ao bombardeamento de electrões que encorpam 
as cescargas parciais, crescem de preferência segundo 
as linhas de campo eléctrico numa estrutura amorfa ou 
pelo menos parcialmente cristalina com as cristalites 
desordenadamente orientadas. Se existir uma orientação 
preponderante das cadeias moleculares sobrepõe-se ao 
efeito do campo eléctrico o efeito de orientação macro- 
molecular [5]. 

Pela análise teórica e experimental da acção do cam- 
po eléctrico no crescimento do arvorejamento verifica-se 
que se a distância disruptiva for muito grande o fenó- 
meno progride num forte adensamento de canais num 
volume limitado por um contorno com formas típicas 
distintas (fig. 6 a): no caso de uma heterogeneidade 
elíptica o arvorejamento é alongado, exactamente como 
as linhas de campo, e num eléctrodo hiperbólico tem 
um aspecio arredondado, em acordo com uma maior 
abertura das respectivas linhas de campo. 

O contorno do arvorejamento, relativo aos extremos 
dos ramos de maior comprimento, é definido pela ener- 
gia dos electrões que atingem esse limite, sob a influên- 
cia da orientação e intensidade do campo eléctrico, da 
ordenação macromolecular do sólido e da pressão dos 
gases resultantes da decomposição química da substân- 
cia polimérica. 


! ! o ou de ds ds 4 


O x [mus] sima) mo 


sjmm] e 
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Fig. 6 — Comparação do arvorejamento de polietileno nas configurações elíptica e hiperbó- 
lica: r= 0,01 mm, s=3 mm, d=4,2mm; linhas de campo (a traço-traço) es- 
timadas a partir das linhas equipotenciais (a traço contínuo fino); linhas isocampo 


(a traço contínuo grosso). 


a) Arvorejamento na configuração elíptica (à direita) e hiperbólica (à esquerda) 


após 15 minutos a 30 KV. 


b) Distribuição do campo (à esquerda) e arvorejamento (à direita) mo modelo 


hiperbolóide-plano. 


c) Arvorejamento (à esquerda) e distribuição do campo (à direita) no modelo 


elipsóide-plano. 
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Mas a forma desse contorno é esclarecida pela con- 
jugação das linhas de campo (que definem a direcção 
do campo) com as linhas isocampo (que determinam a 
sua intensidade). A princípio pensámos que o arvoreja- 
mente deveria atingir exactamente as superfícies isocam- 
po, pois correspondem a pontos com igual intensidade 
do campo eléctrico [6]. Mas observando as linhas de 
campo no caso do eléctrodo elíptico (fig. 6 c) conclui-se 
que as partículas aceleradas são desviadas paralelamen- 
te ao eixo não podendo alcançar os pontos das linhas 
isocampo afastadas do eixo. O andamento das curvas 
isocampo esclarece por seu lado a forma oval do con- 
torno pois para duas trajectórias paralelas ao eixo a in- 
tensidade do campo é tanto menor quanto mais afasta- 
das se encontram as partículas do eixo sobre o mesmo 
plano transversal; a este efeito de redução da intensidade 
do campo sobrepõe-se o efeito de curvatura da trajectória 
(com maior cedência de energia). De facto na base do 
esclarecimento da forma do contorno de arvorejamento 
reside a energia das partículas no momento em que coli- 
dem com as unidades estruturais das macromoléculas, 
isto é, a energia inicial na ponta do eléctrodo (onde o 
campo eléctrico é máximo) adicionada com a energia 
cinética adquirida pelo movimento até à colisão, que é 
tanto mais reduzida quanto maior for a energia cedida 
durante o percurso por colisões com outras partículas 
e devido à curvatura da trajectória por influência do cam- 
po. Assim se compreende que quanto mais afastadas do 
eixo forem as trajectórias das partículas aceleradas me- 
nor é o seu comprimento, e daí o não se atingirem as 
curvas isocampo. 

O nosso pensamento inicial encontra apoio, em cer- 
ta medida, no eléctrodo hiperbólico (fig. 6 b), pois as 
linhas de campo são mais abertas, originando um contor- 
ro circular mais próximo da evolução das corresponden- 
tes linhas isocampo. 


5. CONCLUSÕES 


Observando experimentalmente o arvorejamento de 
polietileno e comparando com a distribuição do campo 
eléctrico em ambas as configurações estudadas consta- 
ta-se um acordo notável com a teoria exposta: 

1. A direcção do campo eléctrico (fig. 4 a) cdeter- 
mina componentes longitudinais e transversais que justi- 
ficam um arvorejamento mais estreito no caso elíptico 
(fig. 6 a). 

2. O andamento das linhas isocampo (fig. 4 b) 
estendendo-se mais no espaço de disrupção, para iguais 
valores da intensidade do campo eléctrico, dão ao eléc- 
trodo elíptico um arvorejamento mais comprido (fig.6 a). 

3. O valor da intensidade máxima do campo eléc- 
trico na ponta do eléctrodo (fig. 4 c), em igualdade de 
condições geométricas (raio de curvatura e distância 
Cisruptiva), sendo muito mais elevada na configuração 
hiperbólica facilita o estudo do arvorejamento com eléc- 
trodos hiperbólicos pois se exigem tensões inferiores 
para se dar início ao fenómeno. 

4. A intensidade do campo eléctrico ao longo do 
eixo do eléctrodo heterogéneo (fig. 4 c) reduz-se em 
breve a valores insignificantes no caso hiperbólico pelo 
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que o arvorejamento se estabiliza num volume limitado 
se a distância disruptiva for suficientemente grande. Este 
limite é uma característica do material e pode ser enca- 
rado como uma grandeza definidora do fenómeno de 
arvorejamento se se normalizar uma configuração geo- 
métrica hiperbólica (raio de curvatura da ponta e distân- 
cia disruptiva bem definidas). No caso elíptico a inten- 
sidade mínima do campo eléctrico é Eg= V/d, pelo 
que se pode dar o progresso do arvorejamento ao longo 
do tempo. 

5. Para um dieléctrico com uma dada espessura 
quanto mais se penetra o eléctrodo elíptico (redução da 
distância disruptiva), maior é a máxima intensidade do 
campo eléctrico na sua ponta; mas esse valor mantém-se 
sempre inferior ao da ponta hiperbólica equivalente (com 
igual raio de curvatura e para a mesma distância disrup- 
tiva) aproximando-se no entanto os respectivos valores 
(fig. 5). Quando porém a distância disruptiva s é redu- 
zida apenas a configuração hiperbólica tem significado 
para o cálculo do campo, tendo-se definido um critério 
de transição através de rectas tangentes às curvas 
log Emax (s) de igual parâmetro r, critério que satisfaz 
absolutamente com a experiência. 

6. Com grandes raios de curvatura a intensidade 
máxima do campo eléctrico é menor que para pequenos 
raios de curvatura em ambos os tipos de configurações 
(fig. 5), mas os valores do modelo hiperbólico mantém- 
"se sempre acima dos correspondentes no sistema elíp- 
tico, o que explica a preferência dos experimentadores 
pelo chamado teste de agulhas, pois com a heterogenei- 
dade hiperbólica é mais fácil desenvolver o arvoreja- 
mento. 
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